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Abstract. Evidence for a disturbed velocity structure in the outer regions of the galaxy NGC 1399 comes from a 
re-analysis of the planetary nebulae data from Arnaboldi et al. (1994). We find a strong rotation along a P.A.=140° 
for R < 140” which is followed by a rapid drop off at larger radii, where the velocity dispersion starts to increase. 
This kinematical behavior can be interpreted as an indication for an interaction scenario. Interaction is advocated 
in previous analysis of the halo regions of this system, based on different dynamical tracers such as globular 
clusters and X-rays, but in all these studies the mass distribution is derived under the equilibrium hypothesis, 
which is not appropriate when an interaction takes place. Here we attempt a non-equilibrium dynamical analysis 
of NGC 1399: with a simple model and under the impulse approximation, we show that the observed kinematics is 
consistent with an energy injection caused by a flyby encounter of NGC 1399 with the nearby system NGC 1404. 
In this approach, we find a mass-to-light ratio, M/Lb=26 Mq/Lq within R = 400”, which is about half of that 
requested when equilibrium is assumed, i.e. M/Lb=56 Mq/Lq. 
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1. Introduction 

Inner parts of galaxy clusters are the birthplaces of 
complex dynamical situations involving interactio ns be- 
tween galaxies. M orpho logical segregatio n (Dre ssier 198C 


Whitmore et al. 1993 , Dressier et al. 1997|) , HI defi- 


light distribution (Ciardullo et al. 1989; Ciardullo et al. 
1991; McMillan et al. 1993[ Ford et al. 1996| ) and share 
the same kinematical behavior where they overlap with 


the integrated stellar light data (Arnaboldi et al. 1994, 
1996 and 1998; Hui et al. 1995). Thus, PNe appear as the 


cienc y of cluster spirals (Haynes et al. 1984 ; Solanes et al. 
2001 ), exten ded e nvelopes in density profiles of cD galaxies 
(Schombert 198C| ), structures in the densi ty and tempera¬ 
ture maps of X-rays (Roettiger et al. |1995|; Irwin & Sarazin 
19961 ; Kikuchi et al. |2000D , are some of the signatures of en¬ 


vironmental effects on cluster galaxies. These interactions 
may also produce disturbed velocity structures in the out¬ 
skirts of cluster galaxies, which is a well established fact 
for spirals. Significant statistics based on extended rota¬ 
tion curves measured by coupling HI with HII emission 
have disclosed a decline at large galactocentric distances 
in contrast with the asymptotic flat behavior for field spi¬ 


rals (Whitmore et al. 1988; Adami et al. 1999; Rubin et 
al. [1999 ) and/or asymmetries (Dale et al. 2001). On the 
contrary, the situation for ellipticals is still unclear owing 
to the lack of classical kinematical tracers in the outer re¬ 
gions of these galaxies. 

A recent development has been offered by the first mea¬ 
surements of the discrete radial velocity fields of globu¬ 
lar clusters (GCs) and planetary nebulae (PNe) in nearby 
early-type galaxies. In particular, PNe follow the stellar 


natural candidates to gauge the stellar kinematics in the 
outskirts of ellipticals. Up to now the limiting factor in 
the use of discrete radial velocity fields has been the small 
number statistics (SNS). In distant galaxies (D>15 Mpc) 
the measured samples typically consist of less than a hun¬ 
dred radial velocities both for PNe and GCs, and some 
doubts may arise on the mass and angular momentum es- 
timate s bas ed on these data alone. Recently Napolitano 
et al. ( [2001 , NAFC hereafter) have shown that, even un¬ 
der SNS regime, discrete tracers do carry some impor¬ 
tant kinematical information and can be successfully used 
for mass and angular momentum estimates and search for 
kinematical signatures of interactions. 

An interesting target in this context is NGC 1399, the 
cD galaxy of the Fornax cluster. Despite its apparent reg¬ 
ular morphology, this object exhibits some peculiarities 
which indicate that its real dynamical status i s stil l un¬ 
clear. The extended stellar envelop (Schombert 1986 ) and 
the unusual overa bundance of globular clusters (Kissler- 
Patig et al. 1999 and references therein), together with 
an increasing velocity dispersion in the outer regions from 
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discrete radial velocity field of both GCs (Grillmair et al. 
1994, Minniti et al. 1998j, K issler-Patig et al. 1999) and 


PNe (Arnaboldi et al. 1994 ), suggest an interaction sce¬ 
nario, which is also supported by X-ray data. The X-ray 
te mpera ture profile of the exten ded gaseous halo (Ikebe et 
al. 1996; Rangajan et al. 1995; Jones et al. 1997) is com¬ 
patible with the velocity dispersion from PNe and GCs. 
Recently, Paolillo et al. (2001) analysed deep ROSAT HRI 
data with an adaptive smoothing technique. They found: 
1) an extended and asymmetric hot gas halo distributed 
out to the cluster scale; 2) a multicomponent density pro¬ 
file which flattens at 1 arcmin with one “shoulder” at 
R = 450", i.e. within the separation between NGC 1399 
and NGC 1404; 3) filamentary structures and holes in 
the intensity map; 4) a displacement of the gas centroid 
with respect the luminous component. All these features 
are qualitatively expected in a tidal interaction scenario 
(D’Ercole et al. 2000| ; Barnes 2000 ). In this context the dy¬ 
namical description of this system is complex : interactions 
are associated to non-equilibrium situations which are dif¬ 
ficult to handle. Previous studies assumed equilibrium to 
infer the mass distribution of NGC 1399 (Bicknell et al. 
1989|; Arnaboldi et al. 1994; Minniti et al. 1998| ; Kissler- 


Patig et al. 199S| ; Saglia et al. j2000 ). They all show an 
increasing mass-to-light ratio ( M/L > 40 Mq/Lq) out¬ 
wards. These values match the local M/L ratios inferred 
for clusters from X-ray measurements, and are consistent 
with those derived from the equilibrium analysis of the 
temperature profile of the X-ray gas around NGC 1399 
(Killeen & Bicknell 1988; Rangajan et al. 1995| ; Jones et 
al. 1997| ). The question is whether these mass estimates are 
realistic: non-equilibrium kinematics can mimic an overall 
highe r mass content and cause mass overestimates (Mihos 
200C). Here we attempt a quantitative non-equilibrium 
analysis to derive the mass distribution of NGC 1399. 
This paper is organised as follow. New estimates of the 
PNe kinematical quantities in the outer halo of NGC 1399 
are obtained in Sect. 2, where some evidence of a dis¬ 
turbed rotational structure is also given. In Sect. 3 the 
PNe rotation curve and velocity dispersion profile are com¬ 
bined with the stellar and GCs kinematics to derive the 
mass distribution via the inversion of the Jeans equa¬ 
tions, under equilibrium hypothesis. We then investigate 
whether NGC 1399 is a system out of equilibrium because 
of ongoing encounters with nearby systems in the core of 
the Fornax Cluster and derive a new mass accordingly. 
Conclusions are presented in Sect. 4. Throughout this pa¬ 
per, we will assume a distance of 17 Mpc for NGC 1399 
(McMillan et al. 1993| ). 


2. Kinematical analysis 

Using Montecarlo simulations in spherical systems, NAFC 
showed that, even under SNS regime, simple fitting func¬ 
tions provide reliable estimates of the position angle, </>zi, 
of the kinematical major-axis, Z\. As in Arnaboldi et al. 
( 1998| ) both a bilinear (BF) and a flat-curve (FC) function 
(Table 1) were interpolated to the radial velocity field, 
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Fig. 1 . PNe radial velocity field of NGC1399. Full dots are 
velocities above the mean velocity, v, open dots are below v. 
The dashed line is the direction of maximum gradient, i.e. the 
kinematical major-axis, derived from the whole sample of 37 
PNe (see discussion in the text). Identification numbers are 
drawn; from Arnaboldi et al. (1994). 


1/rad (X, F ), of the 37 NGC 1399 PNe from Arnaboldi et 
al. ( [1994 ), shown in Fig. [I]. Clearly, the trend of the rota¬ 
tion curve, 1/rot(Zl), is not necessarily matched by these 
functions, but NAFC proved that some information can be 
recovered by plain binning of the velocity data. In prac¬ 
tice, they suggest to average the velocities of all the stars 
falling within a given range of \Z1\ and within a fixed 
distance, dZ 1, from the Z1 axis; the bins are chosen in 
order to secure subsamples with 10 or more data points. 
As shown in NAFC this implies relative errors less than 
30% for both the velocity dispersion and rotation veloc¬ 
ity. For NGC 1399 we have chosen three bins in \Z1\ and 
dZl = 100". The binned radial velocities are shown in 
Fig. I (left panel) together with the BF and FC projected 
along Z1 and V TO t in each bin. 

Rotation - The binned V ro t data (Table ||) suggest that 
the PNe system rotates quite strongly inside R « 140", 
while it does not in the outer regions. This behavior was 
not seen by Arnaboldi et al. (1994), who discussed only 
the results from the bilinear fit. Grillmair et al. (1994) and 
Minniti et al. (1998) do not report any rotation of the GCs 
system out to ~ 8', while Kissler-Patig et al. ( |1999 ) found 
evidence of rotation only at large radii (I/ max =153±93 
km s -1 for R > 5' along a P.A.=120°±40°). They at¬ 
tribute this feature to a tidal interaction or merging of 
cluster galaxies passing throughout the core of the Fornax 
Cluster, rather than to the sign of a relative high local an¬ 
gular momentum. We shall discuss this possibility in more 
details later. In any case, the lack of an inner rotation of 
the GC system supports its decoupling from the stellar 
component. 
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Fig. 2. Left panel: radial velocities (open circles), which are selected along a 100" wide strip about the Z1 axis, are plotted 
against distance along this axis. Full c ircles are the mea n velo cities in each bin, while radial velocities from integrated light data 


(D’Onofrio et al. 1995; Graham et al. 199q; Saglia et al. 200C) are indicated with open symbols (triangles, squares, pentagons). 


The solid line is the bilinear fit and the dot-dashed line represents the flat curve fit to the data along Zl. Vertical dashed lines 
indicate the bin size adopted for the kinematical estimates. Right panel: radial velocities (with the related PN ID numbers) 
are plotted as a function of the projected distance from the center, R = \JX 1 2 + Y 2 . Vertical dashed lines indicate the width of 
radial annuli adopted in the kinematical analysis. 


Velocity dispersion - In each Zl bin, the projected velocity 
dispersion is the standard deviation of the sample (e.g. no- 
fit procedure, see NAFC), corrected for the measurement 
errors, <r p = \/S D 2 — where a mea =70 km s _1 from 
Arnaboldi et al. ( |l99d ). Outside R ~ 140", where no rota¬ 
tion is found, the velocity dispersion is obtained in radial 
annuli in order to retain all data (Fig. ||, right panel) and 
improve precision. The results are summarized in Table 
|j. The velocity dispersion has a minimum at R ~ 140", 
where V ro t drops, and increases outwards (see also Fig. ||). 
PN dp values are consistent with the integrated light’s at 
small R and with GCs’ at large R. 

This peculiar kinematical behavior possibly calls for a 
complex dynamical situation. In particular, the evidence 
that decoupled populations, such as PNe and GCs, exhibit 
the same “heating’’ in the outer regions appears peculiar. 
In the following we shall consider different scenarios for the 
dynamical status of the system, to be extensively analysed 
in the next section: 

Case a). The system, at equilibrium, has two decou¬ 
pled dynamical components: an inner, rotating, cold com¬ 
ponent (possibly a disk) and an outer, non-rotating, hot 
component, and the transition region occurs at R ~ 140". 
In principle, we can expect some signatures of this mul¬ 
ticomponent structure in the surface brightness profile. 
Improved photometric studies and a more detailed kine¬ 
matics are needed for any significant treatments of this 
picture. The possibility of an inner extended disk as cold 
component seems, however, unlikely because such large 


disk sizes are not observed in ellipticals (Nieto et al. 1991 


Case b). In a more conservative way, if PNe are a 
system at equilibrium in the gravitational potential of 
NGC 1399, the rise of cr p (i.e. of the pressure gradient) is 
needed to compensate the lack of rotation at large radii. In 
this and the above picture, the external regions are fading 
into the cluster as we can infer from both the luminosity 
and velocity dispersion profiles which are consistent with 
the c luster profiles (Gril lmair et al. 1994 , Arnaboldi et al. 
1994, Kissler-Patig et al. 1999| ). In the following dynamical 
analysis, we will refer to this as Model I . 

Case c). An alternative picture is that the regions at 
radii R > 140" experienced some forms of interaction with 
the cluster. This possibility has been invoked in previous 
works on PNe and GCs kinematics and also in X-ray stud¬ 
ies (Rangajan et al. 1995; Ikebe et al. 199C| : Jones et al. 
1997; Paolillo et al. 2001). In this case there might be ob¬ 


jects which are no more bound to NGC 1399. There are 
three outliers in the distribution of 14ad for the last radial 
bin (Fig. |j). Their large velocities are inconsistent with the 
Gaussian-like behavior expected for the radial velocity dis¬ 
tributions based on small statistical samples[]. Moreover 
they are responsible for the number overdensity in the 
outer regions of the galaxy (see Fig. j||). Were NGC 1399 
a relaxed system, these objects should not be included 
when computing the mass at equilibrium ( Model II . here¬ 
after). The global kinematical quantities computed after 
discarding PN ID=1, 13 and 36 (Fig. B), are listed in Table 


van den Bosch et al. [1994 Scorza et al. |1998|) or exp ected 
from multicomponent models (van den Bosch 199ij ): this 
scenario will not be included in the analysis which follows. 


1 The observed distribution of radial velocity profiles in early 
type galaxies have typical deviati on of 10% from a Gaussian- 


like behavior (Winsall & Freeman |1993| ; Bender et al. |199J). In 
order to appreciate these deviations, m ore t han few hundreds 
of radial velocities are needed (Merritt 1997). 
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Table 1. Global kinematical quantities from the PNe radial velocity field. The first column gives the number of PNe 
considered (see discussion in the text). I4 ys and Vm ax are in km s -1 and gradV are in km s -1 arcsec -1 . 


PNe global kinematical quantities 



Bilinear Fit 

Flat-Curve Fit 

Sample 
37 PNe 
34 PNe 

Vsys 

1502±60 

1521±49 

‘f’zi 

143°±15° 

134°±12° 

gradV 

1.1±0.9 

0.52±0.11 

4>zi 

140°±15° 
137° ±13° 

fmax 

177±91 

107±23 




Fig. 3. Left panel: Distribution of the radial velocities in the last radial bin, Ru s t = 225” (see Fig.^j). Three PNe (ID=1,13,36 
in Fig. jij) are inconsistent with the Gaussian-like behavior of the remaining objects. Right panel: PNe radial number density 
profile. PNe follow the surface brightness profile out to 200”: at larger radii there is a possible excess of PNe. 


Table 2. Projected kinematical quantities (given in km s - 1 ). Errors on the rotation velocity and velocity dispersion 
are SD/VN and SD/[y/2N(l — cr^ lea /S'I? 2 ) 1 / 2 ] respectively. The latter is obtained by propagating the measurement 
errors. The code in the last column is: NF=no-fit procedure, BF=bilinear-fit procedure, R=radial annuli, Z1=Z1 bins. 
For an extensive description of the procedures, see Napolitano et al. ( 2001 ). W Standard deviation of the sample of 
the residual field AV = Vrad — BF (see NAFC). 


PNe projected kinematical quantities 


37 PNe 

34 PNe 


Distance 

Trot 

<r p 

Distance 

Trot 

CTp 

Procedure 

39” 

16±121 

297±108 

28” 

-3±156 

304±110 

NF(Z1) 

101" 

251±77 

244±57 

101” 

225±75 

262±58 

NF(Z1) 

116” 


260±70 

116” 


257±60 

BF(R) W 

156” 


334±77 

157” 


188±47 

NF(R) 

182” 

8±107 

360±90 

173" 

-80±70 

195±49 

NF(Z1) 

225" 


460±96 

224” 


150±39 

NF(R) 


|l|, and the rotation velocity and velocity dispersion values 
in Table f|. 

Case d). Finally, we assume that the outer halo of 
NGC 1399 is not at equilibrium ( Model III , hereafter). 
This picture is treated separately in the next Section. 


with larger samples (Hui et al. 1995), one solves the Jeans 
equations in the equatorial plane and obtains the mass 
distribution in some analytical form. NGC1399 is a nearly 
round (E0-E1) galaxy. Its dynamics depends on the ellip- 
ticity of the potential, which is much rounder than its den- 


3. Dynamical analysis 

By combining PN and GC data with the inner stellar kine¬ 
matics it is possible to infer the mass distribution out 
to the halo regions of NGC 1399^. As in previous work 

2 This is strictly true for PNe because they share the same 
dynamics than the inner stellar population. As stressed by 


Saglia et al. (200C), GCs number density profile as function 
of radius follow a powe r law with a slope which is 0.28dex 
shallower (Forbes et al. 1998 ). Therefore, we expect GCs to 
have a velocity dispersion « 1.08 larger than stars and PNe; 
when they are used together with the stellar kinematics the 
velocity dispersion has to be corrected for this factor. 
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sity distribution (e$ ~ e p /3; Binney & Tremaine 1987). 


Therefore, spherical symmetry is adopted for this system. 

For sake of simplicity, we assume isotropy for the ve¬ 
locity ellipsoid in the Jeans equations, accordingly with 
previous works based on discrete radial velocity fields in 
different early-type systems (Grillmair et al. 1994, Hui 
et al. 1995, Arnaboldi et al. 


Kissler-Patig et al. 
1999). Indeed, isotropy is not excluded by previous dy¬ 


namical studies, based on integrated light data: i) Bicknell 
et al. ( 1989| ) found a variable anisotropy parameter within 
100" around a central null value for a constant mass-to- 


light ratio, M/Lb = 14; ii) Saglia et al. (200C) found 


that a wide range of anisotropy parameter, including the 
isotropic case, are compatible with the dispersion profiles 
from PNe and GCs. On the equatorial plane (2 = 0), the 
isotropic axisymmetric Jeans equation 


dp(R,0)crR(R,0) p(R, 0) 


dR 


R 


t (R) = -p(R, 0) 


<9<L(R, 0) 

OR 


( 1 ) 


is equivalent to the Isotropic Radial Jeans Equation 
(IRJE) in spherical coordinates. Using p{R , 0) = p(r) from 
the light profile with a fixed M/L, U-ot from the observa¬ 
tions and adopting an analytical form for the potential, 
Eq. (ii;) can be solved and the mass distribution derived 
from the velocity dispersion data. In more details: 

1) under spherical symmetry, we compute p(r) by depro- 
jecting the Killeen & Bicknell’s (1988) surface brightness 
profile 


E b (R) =E 0 ( 1 +J- 


( 2 ) 


with E 0 = 16.73 ± 0.05 mage arcsec -2 , R s = 3.0" ± 0.2" 
and a = 1.72 ± 0.05 and the mass-to-light ratio, M/Lb = 
12 Mq/Lq (we adjust the M/L from the same authors to 
our adopted distance). This profile is consistent with the 
PN e num ber density distribution (see Fig. |j). Saglia et 
al. (2000) performed an accurate deprojection of the light 
distribution and a dynamical modeling for NGC 1399, ob¬ 
taining a M/Lb = 10 Mq/Lq for the luminous compo¬ 
nent inside R = 60", where the influence of the dark halo 
starts. Since they do not provide us with the stellar den¬ 
sity distribution in an analytical form, we prefer to use 
the analytical expr ession from Eq. (||) (as was done by 
Kissler-Patig et al. 1999t for instance). 

2) Due to the lack of detailed data, a simplified model 
of Vrot(-R) is assumed: we have a linear regime out to 
R ft: 140", then V ro t drops to zero. This artificial model is 
deprojected assuming circular orbits and no correction is 
adopted for inclination (the reason is that the peak veloc¬ 
ity is already large enough to accommodate a high incli¬ 
nation) . 

3) We consider a dark mass contribution to the potential 
by a simple pseudo-isothermal sphere: 


Pd = 


2nG(r 2 + r)j) 


(3) 


where <t< j and r<i are free parameters. In the IRJE: 

d4>(r) _ M\ um (r) + M d (r) _ 
dr r 2 

= —T-[ (TiOO + Pd{r'))r , 2 dr' (4) 

r J o 

where 7 = M/Lb and j(r) is the deprojected luminosity 
profile; 

4) Od and are assigned by matching the velocity 
dispersion solution of the IRJE, once projected along the 
line-of-sight, to the observed velocity dispersion estimates 
by minimizing the x 2 - 

This procedure implicitly assumes that the system is at 
equilibrium. 


3.1. The equilibrium analysis: mass distribution 


Model I: We consider the complete sample of 37 PNe plus 


the GCs data from Kissler-Patig et al. (1999). The fit to 
the velocity dispersion profile is shown in Fig. ^ with the 
confidence regions for the best-fit parameters (ad=430 
km s -1 , rd=78"). The mass distribution is shown in Fig. 
[|: it gives M/Lb = 56jli| Mq/Lq within 400", where the 
quoted errors are related to the extreme values of Cd and 
rd in 68% confidence region and produce the dot-dashed 
dispersion models in Fig. 0. 


Model II: Here only 34 PNe (the three PNe (ID = 
1, 13, 36) were excluded) are considered to be at equi¬ 
librium. Fig. [| shows the fit to the velocity dispersion 
estimates with the related confidence regions in param¬ 
eter space. The best-fit parameter are <7d=280 km s -1 
and rd=44". In this case, we introduced a tidal radius, 
i?t=360", where the stellar distribution is truncated, 
which has improved the fit to the data. This possi¬ 
bly indicates the most external radius of NGC 1399, 
while the cluster potential dominates at larger radii. 
The derived mass distribution is also shown in Fig. ^ 
in this case we obtain M/Lb = 3 S/^Mq/Lq within 400". 


Despite the oversimplified approximation for the rota¬ 
tion curve, the above analysis provides the mass distribu¬ 
tion needed to match the velocity dispersion profile in the 
two extreme rotation regimes (i.e. rigid rotation within 
R«140", null rotation outside). The link between the two 
regions calls for a more detailed kinematics at the interme¬ 
diate radii. The mass distribution from our analysis is in 
perfect agreement with estimates from Saglia et al. ( 2000 ) 
who found M=1.2-2.5x 10 12 M S within R=400". 


3.2. The non-equilibrium scenario for NGC 1399 

The peculiar kinematics and the temperature profile of 
the hot gas from the X-ray observations may sugges t that 
NGC 1399 is not a relax ed sy stem. Minniti et al. ( 1998 ) 
and Kissler-Patig et al. ( 1999| ) use the metallicity distri¬ 
bution and specific density of GCs to claim evidence for a 
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Fig. 4. Velocity dispersion profiles observed and modeled and the confidence regions for the model parameters. Upper left box: 
the PNe estimates in Z1 bins (filled pentagons) and radial annuli (filled circles) are plotted with the mean stellar kinematics 
in spatial bins (filled triangles) obtained from different authors (plotted with different symbols) and GCs estimates (open 
pentagons, corrected for the factor 1.08, see discussion in the text). Under equilibrium hypothesis (Model I), the solid line is the 
best-fit (crd=430 km s -1 , rd=78", 98% significance level ) to the data, the dashed line is the velocity dispersion profile derived 
from the temperature of the X-ray gas (Jones et al. 1997). Dash-dotted lines correspond to the extreme parameter values in the 
68% confidence region. Upper right box: we show the confidence regions at 68%, 90%, 95% in the parameter space. Lower left 
box: velocity dispersion estimates for the 34 PNe sample. Solid and dash-dotted lines as above: here they are computed for the 
Model II as best-fit to the data (<7d=280 km s -1 , rd=44", 98% s.L). In this case no GCs data are used. Lower right box: the 
68%, 90%, 95% confidence levels for Model II. 


Table 3. Projected kinematical quantities (units are km 
s -1 ) for the PNe subsample within 140". 


PNe projected kinematical quantities 


Model III 


Distance 

Kot 


Procedure 

33 

109±104 

221±73 

NF(Z1) 

88 


244±63 

BF(R) 

94 

200±72 

192±51 

NF(Z1) 

126 


146±40 

BF(R) 


recent interaction with the nearby companion NGC 1404. 
If such interaction occurred, the most spectacular effect 
might be the strong heating of the intergalactic gas. The 
hot gas in NGC 1399 has a temperature which, once 
converted into a velocity dispersion, seems systematically 
larger by ~30 km s -1 than that of the stellar popula¬ 


tion (see the long-dashed curve in Fig. [J). Moreover, the 
decline of the rotation curve of NGC 1399 is similar to 
those observed in spiral systems and reproduce d by s ingle 
and/or multiple galaxy encounter (Kauffmann 1999; Salo 
& Laurikainen [2000 ). 

We shall then explore a model where NGC 1399 is 
out of equilibrium, and the stars in the outer halo 
are re-arranging themselves after a flyby encounter with 
NGC 1404. NGC 1399 and NGC 1404 have a small pro¬ 
jected distance, b = 9', and a large relative radial velocity, 
V = 522 km s' 1 (V1399 = 1425 ± 5, Graham et al. [L998, 


V1404 = 1947 ± 5, D’Onofrio et al. |1995| ). The simplest ap¬ 
proach is to use the general scheme of the impulse approxi¬ 
mation (Binney & Tremaine 1987), where NGC 1399 is the 
perturbed system and NGC 1404 is the perturber. In this 
approach, the encounter does not modify the potential of 
the perturbed system, and the kinetic energy in the inner 
regions, while the outer regions (R > 140") experienced 
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Fig. 5. The mod eled mass distributions for Model I (left) and Model II (right). Dashed line is the luminous mass using Killeen 
& Bicknell (1988) luminosity profile with M/Lb = 12M@ /L@, dash-dotted line is the dark halo mass from Eq. (^[) where the 
best-fit parameters were adopted, and the solid line is the total mass, Mtot- 


1 ) an energy injection, 2) the disruption of the stream¬ 
ing motions, i.e. of the angular momentum (Sugerman et 
al. [2000 ), both re-distributed in the random motions, and 
3) no mass loss (we will discuss this assumption later). 
Following this scheme, the post-encounter kinetic energy 
of the perturbed system, Ea n , can be written as 


Ea n — E[ n + A E 


(5) 



where E- ln is the kinetic energy before the encounter, that 
is at equilibrium, and A E is the variation in the ki¬ 
netic energy induced by the encounter. In Eq. (|5|), E^ n 
is associated to the observed random motions only, which 
are the result of the energy injection and the disrup¬ 
tion of the streaming motions; it can be computed as 
£)} n = 3/2 Mo' 2 hs . Let us estimate E- m and A E. 

The equilibrium before the encounter - E- m writes as 


E in = 1/2 MV? ot + 3/2 Mal v 


( 6 ) 


where V TO t and er eq are derived from the kinematical data 
of both the integrated light and PNe in the range of 
radii where the system is assumed to be still at equi¬ 
librium (R < 140"). In particular, for the PNe system 


-l 

-l 


$zi = 131° ± 7°, 
(BF). Here we em- 


we obtain I4 ys = 1465 ± 73 km s 
gradV = 2.1 ± 0.3 km s -1 arcsec 
phasize that the 14ys estimate, from the inner PNe sub¬ 
sample, is in good agreement with the stellar light data. 
The estimates in Table show a positive residual in the 
outer regions, i.e. in the same direction of the assumed en¬ 
counter between NGC 1399 and NGC 1404. This effect is 
expe cted in the impulse approximation (Aguilar & White 


1985), and supports our assumption about the inner equi¬ 


librium. 



Fig. 6. Model III. Top panel: rotation velocities from stars 
(open circles) and PNe (full circles) with the modeled rotation 
curve from Eq. (|7(). Bottom panel: velocity dispersion estimates 
and the best-fit to the data within 140". Symbols have the same 
meaning as in Fig. [ij. 


V TO t and o p are listed in Tabic ^ and shown in Fig. [| 
In the left panel V TO t is interpolated with a flat rotation 
curve of the formf^] 


V rot (R) 


VqR 

s/R 2 + Rn 


(7) 
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Fig. 7. Model III. The total mass distribution (heavier solid 
line) for Model III is compared with mass distributions from 
Model I (short-dashed line) and II (dotted line). Heavier long- 
dashed line and dot-dashed have the same meaning than as in 
Fig. 1 The contribution of the dark halo becomes important 
for R > 60" in all models. 


where Vo = 214 km s _i and i?o = 67". This is used in 
the IRJE to fit er p as in Sect. 3 and obtain the mod¬ 
eled <r eq . The “ best-fit” is shown in Fig. ^ (right panel): 
the parameter are aa=230^74 km s _1 , ra = 67t 2 s arcsecs 
(P(x 2 > Xo) quoted errors are 68% c.l.) and the 

mass distribution is in Fig. [f]. The mass-to-light ratio is 
M/Lb = 26^g Mq/Lq within 400", where the errors refer 
to the extreme values of the 68% confidence region. 

The injected energy AE - Following Binney & Tremaine 
( 1987 ), we estimate the injected energy using the tidal ap¬ 
proximation equation 


A E = 


AG' 2 M/Mi — 
3 b A V 2 T 


( 8 ) 


analysed in their work (for our model we obtain (3 =0.8 
and p =4, see r 1 / 4 model in Figs. 2 and 4 in Aguilar & 
White ( 1985 )). The corresponding AM/M is about 4%, 
which is well within our M/L uncertainty: this means that 
our assumption of null mass loss is consistent with our 
model precision. Similar conclusions are found by com¬ 
paring our results with Funato & Makito (1999) : in this 
case the expected A M/M is about 5% (for V/a=3 and 
p/v T =2 for their Hernquist model, Fig. 11 c) and d)). 

Eq. (|8|) is the most direct way to estimate the injected en¬ 
ergy, but not the most accurate. Aguilar & White (1985) 
pointed out that this “distant encounter” approximation 
can overestimate the real amount of energy change for 
a large range of impact parameters. From inspection of 
their Fig. 4, however, we can expect some effect for our 
encounter geometry : [3 =0.8 and p =4 imply an overes¬ 
timate of ~ 0.02 in A E/E. We will consider this source 
of uncertainty in the following discussion about the error 
budget of this analysis. 

In conclusion, we stress that, within the flyby encounter 
scenario, we have neglected possible effects of the rela¬ 
tive nearby system, NGC 1380. This galaxy has a pro¬ 
jected distance of 38', an equally large relative velocity, 
¥'=452 km s _1 , with respect NGC 1399 and is as bright 
as NGC1404 and possibly as massive. This large projected 
distance, anyway, would produce a negligible AE' in a 
present tangential encounter (a flyby encounter in the past 
cannot be excluded, as suggested by Kissler-Patig et al. 


1999). 


3.2.1. Recovering <r 0 b s : consistency with the 
non-equilibrium scenario 

We aim at verifying whether the velocity dispersion profile 
in the outer regions of NGC 1399 can be explained in terms 
of Eq. (|^): we will do it for the last Z1 bin, Z\ ast = 182". 
Using the value for AE in the definition of the post¬ 
encounter energy, Eg n we obtain 


where M\ and M 2 are the total masses of NGC 1399 and 
NGC 1404 respectively, b is the impact parameter, V is 
the relative velocity, r 2 is the squared mean radius of the 
perturbed system and G is the gravitational constant. We 
initially assume a tangential encounter, where b = 9' is the 
projected separation between NGC 1399 and NGC 1404, 
and the relative velocity is then given by the difference of 
the object systemic velocities, V = 522 km s _1 . In Eq. (||), 
we adopt M\ = 1.7 x 1O 12 M 0 and (r 2 ) 1 / 2 = 40 kpc from 
our mass model. M 2 is computed from NGC 1 404 t otal 
luminosity Lb = 1-84 x 1O 1O L 0 (de Vaucouleurs 1991 and 
D=17 Mpc) and a constant M/Lb=17 Mq/Lq, similar 
to NGC 1399. Finally we obtain AE = 9 x 1O 15 M 0 km 2 
s -2 , which is ~7% of the total energy. 


Following Aguilar & White (1985), this is consistent with 
the energy change expected in similar flyby encounters 


different ellipticals (NG C 5128: Hui et al.__1995; GC 3115: 
Capaccioli et al. 


Arnaboldi et al. 1993). 


1993 ; Emsellem et al. 


1999 


NGC1316: 


lMo 2 ohs = ^MV r l t + ^Ma 2 eq + AE. (9) 

M is the mass content of the region where the injected 
energy is redistributed which, in our case, is the last bin 
along the Z1 axis. Solving with respect to <r 2 bs , we have 


2 _ _ t/2 

°obs ^ 'rot 


9 2 AE 

(7 H-. 

eq 3 M 


( 10 ) 


Using the values of Vrot and er e q extrapolated from the 


equilibrium configuration before the encounter, Eq. (10) 
gives dobs = 264 ±43 km s _1 , which is consistent with the 
estimates from the kinematical data for both PNe and 
GCs (within the errors at Zi as t)- Here the quoted errors 
are obtained by propagating the 1 <t model uncertainties 
on a eq (~ 35 km s _1 ), M(r) (#M/M=0.3) and V rot (~ 
10 km s~ 4 ) at Ziast considered as dependent quantities. 
This result is consistent, within the estimated error, with 
the possible overestimate introduced using Eq. (0). In this 
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Fig. 8. cr 0 bs as a function of b and V: contours are from 270 
km s' 1 to 450 km s -1 with a step of 20 km s _1 , the cross is the 
location of the crtan estimate. See the discussion in the text. 


sense our A E estimate can be considered reasonably ac¬ 
curate. 

Larger uncertainties come from our assumption on the en¬ 
counter geometry. The tangential configuration is neither 
the most obvious, nor the most favorable: it is not diffi¬ 
cult to envisage configurations where the relative velocity 
between the systems is lower than their relative velocity 
along the line of sight and/or their projected distance is 
larger than their distance at closest approach. In this lat¬ 
ter case the amount of injected energy increases together 
withcTobs- In Fig. ||, er 0 b s from Eq. (10) is shown as a func¬ 
tion of b and V: the contours cover the ±1RMS range of 
the CTi as t in the last bin, the cross is our tangential es¬ 
timate, cJtan? and the grey region is the overlap between 
FIRMS values of both cri ast and a tan - The tangential esti¬ 
mate should be considered as an indicative result, because 
a wide range of configurations can match the observed dis¬ 
persion. Conclusive results can be obtained by acquiring 
information on the geometry of the encounter. This uncer¬ 
tainty can be solved with more extended surveys of PNe 
in the NGC 1399-NGC 1404 system which can possibly 
disclose the orbit of an ongoing encounter (Moore et al. 
19961 ). 


Finally, if our assumptions about the inner equilibrium 
are correct, this model has two consequences: 1) the mass 
distribution given by the equilibrium before the interac¬ 
tion (Fig. [?]) is a reasonable description of the actual mass 
content in NGC 1399, since the high-speed encounter does 
not imply a mass modification; 2) the observed kinemat¬ 
ics can be explained in term of an energy injection rather 
than a pure pressure gradient needed to balance a steeply 
increase in the mass-to-light ratio. The non-equilibrium 
hypothesis leads to an outer M/L ratio for NGC 1399 
which is consistent with the typical values for giant E 
galaxies obtained using the stellar (including PNe) and 
GCs kinematics when compared at the same galactocen- 


tric distances (NGC5128: Harris et al. 1988; Hui et al. 


2001 


1995| ; NGC1316: Arnaboldi et al. 1998| Goudfrooij et al. 
NGC3115: Emsellem et al. 199S; see also Magorrian 


& Ballantyne 2001 for a collection of 18 early-type galax¬ 


ies . 


4. Discussion and Conclusions 

So far the discrete radial velocity fields of PNe and GCs 
are the only viable diagnostics to extend the dynami¬ 
cal studies into the outer regions of spheroidal galaxies. 
Furthermore PNe, being associated to the stellar popula¬ 
tion, allow us to extend the stellar kinematics out to the 
halo regions. Despite the small number statistics, PNe and 
GCs can be used efficiently for “normal” systems, while 
they suffer some limitations for the “disturbed” systems, 
where a more detailed kinematical information is needed. 
In this latter case, we can still hope to derive some indica¬ 
tions of the actual dynamical state of their parent galaxies. 
This is possibly the case of NGC 1399, the cD galaxy in 
the center of the Fornax cluster. Here, the re-analysis of 
the PNe data from Arnaboldi et al. (1994), by mean of 
the procedure tested in Napolitano et al. ( [2001 ), shows a 
peculiar rotational structure: V TO t has a peak of ~ 250 km 
s -1 at R ~ 140" and then declines. At the same distance 
from the center, the velocity dispersion has a minimum: 
it decreases out to R = 140" and then increases. This 
could be a kinematical signature of dynamical substruc¬ 
tures (a truncated disk, for example) in the inner regions 
which is not highlighted in published photometrical stud¬ 
ies. Anyway, the fact that GCs show the same dispersion 
profile as the PNe, even if they are a decoupled population 
(they do not rotate in the regions where they overlap with 
the PNe data), make this multicomponent scenario inap¬ 
propriate to explain the observed “heating” of the outer 
kinematics. This heating possibly calls for a different dy¬ 
namical justification. 

In a conventional way, assuming equilibrium, we derived 
the mass distribution of NGC 1399 accounting for the ob¬ 
served kinematics and found a large dark-mass content 
with M/L=56 (Model I) at R=400", fully consistent with 
previous analyses. In a different approach, the peculiar 
kinematics is interpreted as the signature of an interaction 
situation. Such an interaction scenario has been invoked in 
previous work to account for the peculiarities of this sys¬ 
tem: the extended stellar envelope (Schombert 1986), the 
unusual overabundance of globular clusters (Kissler-Patig 
et al. 1999| and references therein), and their metallicity 
distribution (Minniti et al. 1998| ) . The interaction scenario 
is also compatible with X-ray observations: Paolillo et al. 
(2001) have recently analysed new data from ROSAT HRI 
deep observations and found features (holes, filamentary, 
shoulders in the density profile) which are qualitatively 
expected in a tidal interactio n scen ario as shown by ana¬ 
lytical works (D’ Ercole et al. |2000| ) and by N-body sim¬ 
ulations (Barnes 200C). These independent evidences of 
ongoing interaction make any equilibrium analyses in the 
outer regions of this galaxy quite uncertain: in order to be 
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conservative, this is done, however, in Model II, where we 
find M/L=33 at R=400". 

Here we have presented a simplified model as an attempt 
to approach non-equilibrium situations in a quantitative 
way. The lack of a detailed kinematics implies assumptions 
on the geometry of the encounter which affect the final re¬ 
sult and therefore the dynamical situation of NGC 1399 
requires further investigations. However, this simple dy¬ 
namical analysis gives two interesting indications: 1) the 
observed kinematics is coherent with a scenario of non¬ 
equilibrium due to a galaxy-galaxy flyby encounter; 2) by 
considering non-equilibrium, the total mass content re¬ 
quired within R=400" is ~30% less than for the equilib¬ 
rium model (Model I) and implies M/L=26. 
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